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We summarize the main features of the spectroscopy, production and decays of the J PC = 1" 
tetraquarks in the bb sector, concentrating on the lowest state called (10890). The tetraquark 
framework is used to analyze the BaBar data on the e + e~ — > bb cross section (Rf, energy 
scan) between y/s = 10.54 and 11.20 GeV and the Belle data on the processes e + e~ — > 
Y(1S)tz + tz~ ,Y(2S)n + 7Z~ near the peak of the T(5S) resonance. The BaBar R), energy scan is 
consistent with an additional state at a mass of 10.90 GeV and a width of about 28 MeV, in broad 
agreement with the state Ti(10890) GeV seen by Belle in the exclusive final states. We argue 
that the decay widths and the dipion invariant mass distributions measured by Belle are naturally 
explained by the tetraquark interpretation of Fj, (10890). 
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1. Introduction 

Experiments at the B factories and Tevatron have in the past several years revived the in- 
terest in the spectroscopy of the Quarkonium-like exotic states. Labeled tentatively as X, Y and 
Z, due to a lack of consensus on their interpretation, they have masses above the open charm 
(DD) threshold, with the X(3872) being the lightest and y(4660) the heaviest state observed so 
far [1], There is also evidence for an exotic ss bound state Fy(2175) having the quantum numbers 
J PC = 1 , first observed in the initial state radiation (ISR) process e + e~ — > 7isr/o (980) 0(1020), 
where /o(980) is the ++ scalar state. In the bb sector, Belle [2] has observed enhanced produc- 
tion for the processes e + e~ — > Y(IS)ti + 71~ ,Y(2S)ti + 71~ ,Y(3S)ti + 71~ in the e + e~ center-of-mass 
energy between 10.83 GeV and 1 1.02 GeV, which does not agree with the conventional Y(5S) line 
shape [3]. The enigmatic features of the Belle data are the anomalously large decay widths for 
the mentioned final states and the dipion invariant mass distributions, which are strikingly differ- 
ent from the conventional QCD expectations for such dipionic transitions. A fit of the Belle data, 
using a Breit-Wigner resonance, yields a mass of 10888+^(stat) ± 1.2(syst) MeV and a width of 
30.7+^ (stat) ±3.1 (syst) MeV [2]. This particle is given the tentative name Y h (10890). In [4, 5], 
Yb( 10890) is interpreted as a bb tetraquark state, which is a linear superposition of the 7 PC = 1 
flavour eigenstates Y^ = [bd][bd] and Y^ bu ] = [bu][bu]. The mass eigenstates Yn,n (for the lighter) 
and Y\ b u (for the heavier) of the two are almost degenerate, with their small mass difference aris- 
ing from isospin-breaking [6]. A dynamical model for the decay mechanisms of Yy{ 10890) and 
the final state distributions measured by Belle was developed in [4] and refined in [5], yielding 
good fits of the Belle data. One anticipates that 7/^(10890) is also visible in the energy scan of the 
e + e~ — > bb cross section, which was undertaken by the BaBar collaboration between yfs = 10.54 
GeV and 1 1.20 GeV [7]. A fit of the BaBar data on /^-scan is consistent with a structure around 
Yj,{ 10890) and yields a better # 2 /d.o.f. than the fits without the tetraquark states. More data are 
required to resolve this and related structures in the Rj, line shape. This contribution summarizes 
the work done in [4, 5, 6] interpreting the Belle [2] and BaBar [7] data in terms of the bb tetraquark 
states. 

2. Spectrum of bottom diquark-antidiquark states 

The mass spectrum of tetraquarks [bq] [bq 1 ] with q = u, d, s and c can be described in terms of 
the constituent diquark masses, iuq = miu, spin-spin interactions inside the single diquark, spin- 
spin interaction between quark and antiquark belonging to two diquarks, spin-orbit, and purely 
orbital term [8], i.e., with a Hamiltonian 

H = 2m Q + H { s f ] +Hf s Q) +H SL +H LL , (2.1) 

where: 

=2(^) 3 [(S*.S ? ) + (Sg.S 3 )], 
hJs^ = 2(^fc^)(Sfo • + • S q ) +2J^ h j J (Sb ■ Sg) + 2.y(f cl q{S q ■ S^), 

H SL = 2A e (S Q -L + S s -L), H LL = B Q LQQ{L ^ a + l) . (2.2) 

Here (-#2,^)3 is the coupling of the spin-spin interaction between the quarks inside the diquarks, 
J^l,q are the spin-spin couplings ranging outside the diquark shells, Aq is the spin-orbit coupling 
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of diquark and Bq characterizes the contribution of the total angular momentum of the diquark- 
antidiquark system to its mass. 

The parameters involved in the above Hamiltonian (2.2) can be obtained from the known me- 
son and baryon masses by resorting to the constituent quark model [9]: H = Y,i m i + Hi<j2J%j{Si • 
Sj), where the sum runs over the hadron constituents. The coefficient Jfij depends on the flavour 
of the constituents i, j and on the particular colour state of the pair. Using the entries in the 
PDG [3] for hadron masses along with the assumption that the spin-spin interactions are inde- 
pendent of whether the quarks belong to a meson or a diquark, the results for the masses corre- 
sponding to the tetraquarks (q = u,d,s,c) were calculated in [6]. The lowest eight 1 
tetraquark states [bq] [bq] (q = u,d), which are all orbital excitations with Lqq = 1, have the follow- 
ing spin and orbital angular momentum eigenvalues: Y^l (Sq = 0, Sq = 0, Sqq = 0, Lqq = l), 

Y [bc,} ( s Q = l > S Q = °> S QQ = h Lqq = l), Y* q] (Sq = 1, S Q = 1, S q q = 0, L q q = l), and 

Y [bq\ ( S Q = X ' S Q = X ' S QQ = 2 > L QQ = 0- Identif y m g the lowest lying J PC = 1 state F^] with 
the Yb( 10890) measured by Belle, and using the estimates for the other parameters entering in 
Eq. (2.2), fixes the diquark mass iuq = mn, q ] =5.251 GeV. The uncertainties on the masses of the 

In) 

other six states M Y in = 2,3,4) are higher, as they depend in addition on the mass-splittings be- 
tween the good and bad diquarks, A = t}1q{Sq = 1) — iiiq(Sq = 0), estimated as A ~ 200 MeV [10, 
11]. The central values of their masses are: M { y ] =11133 MeV, = 1 1257 MeV, and = 

1 1227 MeV. Assuming isospin symmetry, the states Y^\ and Y^l are degenerate for each n. Includ- 
ing isospin-symmetry breaking lifts this degeneracy with the mass difference between the lighter 
and the heavier of the two states estimated as M(Y^V) -M(Y^l) = (7 ±3)cos(20) MeV, where 

6 is a mixing angle and the mass eigenstates are defined as: Y^l = cosQFrS + sin OY^j, and 
Y [bh] = ~~ sm ^[i'rf] ^ cos ® Y {bu}[^- ^ ne resu l tm g mass differences are small. However, depending 
on 6, the electromagnetic couplings of the Y^l and Y^h ma ^ tum out *° ^ e significantly different 
from each other, and hence also their contributions to Ry. 

3. Decay Widths of ^(10890) and other J PC = 1 tetraquarks 

As the masses of all the eight J PC = 1 [bq] [bq] tetraquark states lie above the thresholds for 

(n) (*) -(*) 

the decays F^j — > B q B q , they decay readily into these final states. For the n = 3 state (having a 

(3) 

mass of 1 1257 MeV), also the decay Yh\ — > A^A/, is energetically allowed. In [6], the decay widths 

r(7^ — > B q *^B q *^) have been estimated (up to a tetraquark hadronic size parameter k) in terms of 

the corresponding partial decay widths T(T(5S) — > B q *^B q *^), which can be calculated with the help 
of the entries in the PDG [3]. Specifically, the following relations are assumed 

K 2 {B + B-\H\Y^ ] u] ) = fc 2 (B°B° \H \ ) = (fl + fl-|#|T(5S)> = (fl°fl°|#|T(5S)) , (3.1) 

and likewise for the BS* and B*S* decays. Noting that the decays (B+B \A\yH), (B°B°\ \H\yW ) 
as well as the decays M^Ud — >■ B s B s ) are Zweig-forbidden, one expects, concentrating on the 
lowest mass state, T(Y^l) ~ 0.4r(T(5S). Using the PDG value [3] T(T(5S) = 110 MeV, we 
get r(F [ | ) 1 jj) = r(Y^\ = (44 ± S)k 2 MeV for the total decay widths. Equating this decay width 
to the measured value of the total decay width r[7 fe (10890)] = 30.7^ (stat) ± 3.1 (syst) MeV 



3 



Tetraquark interpretation 0/3^(10890) and Ry energy scan 



Ahmed Ali 



by Belle [2], one gets K = ^/ff^f = 0.8 ±0.1. This suggests that the tetraquarks F^ have a 
hadronic size of the same order as that of the T(55). The hadronic widths of the other J PC = 1 
tetraquarks are estimated as [6]: r(y^) = 80 ± 16 MeV, r(F^) = 114 ±22 MeV and r(F^) = 
102 ±20 MeV. 

To calculate the production cross sections, we have derived the corresponding Van Royen- 
Weisskopf formula for the leptonic decay widths of the tetraquark states made up of point-like 
diquarks [5]: 



24a 2 |<2 [6u/M] 



m Y b 



(3.2) 



where a is the fine-structure constant, <2[fa<] = ±1/3, Q\bd] = —2/3 are the diquark charges in units 
of the proton electric charge, and \r[^ (0) | 2 = 2.067 GeV 5 [12] is the square of the derivative of the 
radial wave function for Xb(lP) taken at the origin. Hence, the leptonic widths of the tetraquark 
states are estimated as [5] 

r(Y [hd] -> e+e-) = 4r(F [H -> e+e~) « 83 K 2 eV , (3.3) 

which are substantially smaller than the leptonic width of the T(55) [3]. The electronic widths of 
the mass eigenstates Yn,n and Y^m depend, in addition, on the mixing angle 8. 

4. Analysis of the BaBar data on 7^-scan 

BaBar has reported the e + e~ — > bb cross section measured in a dedicated energy scan in the 
range 10.54 GeV and 11.20 GeV taken in steps of 5 MeV [7]. Their measurements are shown 
in Fig. 1 (left-hand frame) together with the result of the BaBar fit which contains the following 
ingredients: A flat component representing the to-continuum states not interfering with resonant 
decays, called A nr , added incoherently to a second flat component, called A,-, interfering with two 
relativistic Breit-Wigner resonances, having the amplitudes A io860>^ 11020 and strong phases, 010860 
and 0no2O> respectively. Thus, 

o(e + e-->bb) = |A, ir | 2 ±|A,.±A 10 860^ 1086 W(M 1 o860,r 1 o86o) 

±Aiio20^ 11020 5W(M 110 2o,r 11 o2o)| 2 , (4.1) 



with BW(M,F) = l/[(s — M 2 ) ± iMT]. The results summarized in their Table II for the masses 
and widths of the T(55) and T(65) differ substantially from the corresponding PDG values [3], 
in particular, for the widths, which are found to be 43 ±4 MeV for the T( 10860), as against the 
PDG value of 1 10 ± 13 MeV, and 37 ± 2 MeV for the T(l 1020), as compared to 79 ± 16 MeV in 
PDG. As the systematic errors from the various thresholds are not taken into account, this mismatch 
needs further study. The fit shown in Fig. 1 (left-hand frame) is not particularly impressive having 
a # 2 /d.o.f. of approximately 2. 

The BaBar /?^-data is refitted in [6] by modifying the model in Eq. (4. 1) by taking into account 
two additional resonances, corresponding to the masses and widths of Y[bj\ an d ^[b,h] ■ Thus, formula 
(4.1) is extended by two more terms 

A Y[b ^MBW(M YM ,r Y[bJ] ) and Ay^/HMBWiMy^Ty^), (4.2) 

which interfere with the resonant amplitude A r and the two resonant amplitudes for T(55) and 
T(65) shown in Eq. (4.1). Using the same non-resonant amplitude A nr and A r as in the BaBar 
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Figure 1: Measured R/, as a function of y/s with the result of the fit with 2 Breit-Wigners described in Fig. 1 
of B. Aubert et al. [BaBar Collaboration] [7] (left-hand frame). The result of the fit of the Ry data with 4 
Breit-Wigners [6] (right-hand frame). Location of the T(5S), T(6S), the tetraquark states (labelled as 

yW) and (labelled as Y^) are indicated. The shaded bands around the mass ofyW andy< 2 ) reflect the 
theoretical uncertainty in their masses. (From [6]). 



analysis [7]. the resulting fit is shown in Fig. 1 (right-hand frame). Values of the best-fit param- 
eters yield the masses of the Y(5>S) and Y(6>S) and their respective full widths which are almost 
identical to the values obtained by BaBar [7]. However, quite strikingly, a third resonance is seen 
in the 7?£-line-shape at a mass of 10.90 GeV, tantalisingly close to the F/,(10890)-mass in the Belle 
measurement of the cross section for e + e~ — > 7^(10890) — > T(IS,2S) 7l + K~, and a width of about 

(2) 

28 MeV In the region around 11.15 GeV, where the Yh{ states are expected, our fits of the BaBar 

(2) 

7^-scan do not show a resonant structure due to the larger decay widths of the states KM. The 
resulting % 2 /&.o.f. = 88/67 with the 4 Breit-Wigners shown in Fig. 1 (right frame) is better than 
that of the BaBar fit [7]. 

The quantity M ee {Y b ) = ~^ee(Y[b,i}) f^ee(Y[b ,h]) 1S given by the ratio of the two amplitudes Ay [bl] 
and Ay [b A| , which also fixes the mixing angle 6 . From the fit shown in the right-hand frame in Fig. 1 , 
one obtains 

@ ee (Y b ) = 1.07 ±0.05, (4.3) 

yielding 

= -19±l° and AM = 5.6 ±2.8 MeV, (4.4) 
for the mixing angle and the mass difference between the eigenstates, respectively. For the mass 
eigenstates and Y^ h h j, the electronic widths T ee (Y[ b jy) and r ee (Y^ bh }) are given by [5] T ee (Q) = 
0.2 K 2 Q(d) 2 keV. With the above determination of fc and 6, we get 

T ee {Y [bA ) = 0.033 ± 0.006 keV and T ee (Y [bM ) = 0.031 ±0.006 keV. (4.5) 



5. Analysis of the Belle data on e + e~ -+Y b -¥ (T(1S),T(25))^ + ^^ 

With the / PC = 1 for both Y b and T(nS), the dipionic final state is allowed to have the quan- 
tum numbers ++ and 2 ++ . There are only three low-lying resonances in the PDG which can con- 
tribute as intermediate states, namely, the two ++ states, /o(600) and /o(980), which we take as 
the lowest tetraquark states, and the 2 ++ qq meson state /2(1270). All three states contribute for the 
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dr/dm„ dT/dco$0 [GeV] dr/dm„ dT/dcos6 [GeV] 




Figure 2: Left-hand frames: Fit results of the Mjj+j,.- distribution and the cos 9 distribution for e + e~ Y{, — > 
Y(lS)n + n~ , normalized by the measured cross section by Belle [2]. Right-hand frames: The same distri- 
butions for e + e~ — > Y(2S)n + n~ . In all figures, the histograms represent the fit results based on tetraquarks, 
while the crosses are the Belle data [2]. The solid curves in the figures show purely continuum contributions. 
(From [4].) 

final state T(lS)n + n^ . However, kinematics allows only the /o(600) in the final state Y(2S)7r + 7r~. 
In addition, a non-resonant contribution with a significant D-wave fraction is needed by the data on 
these final states. This model accounts well the shape of the measured distributions, as shown in 
Fig. 2 for e+e- ->Y h -> Y(lS)n + n- (left-hand frames) and for <? + e~ -> Y b -> r(2S)7T + 7!r (right- 
hand frames). As the decays Yy — > (Y(lS),Y(2S))n + 7i~ are Zweig-allowed, one expects larger 
decay widths for these transitions, typically of 0(1) MeV [5], than the decay widths for the con- 
ventional dipionic transitions, such as T(45) — > Y(IS)ti + 71~ , which are of order 1 keV [3], Further 
tests of the tetraquark hypothesis involving the processes e + e~ — > Yb — > Y(IS)(K + K~ , T]7r°) are 
presented in [5]. 
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